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yA//T7S. We analyze coronagraph observations of a polar jet observed by the Sun Earth Connection Coronal and Heliospheric 
Investigation (SECCHI) instrument suite onboard the Solar TErrestrial RElations Observatory (STEREO) spacecraft. 
Methods. In our analysis we compare the brightness distribution of the jet in white-light coronagraph images with a dedicated kinetic 
particle model. We obtain a consistent estimate of the time that the jet was launched from the solar surface and an approximate initial 
velocity distribution in the jet source. The method also allows us to check the consistency of the kinetic model. In this first application, 
we consider only gravity as the dominant force on the jet particles along the magnetic field. 

Results. We find that the kinetic model explains the observed brightness evolution well. The derived initiation time is consistent with 
the jet observations by the EUVI telescope at various wavelengths. The initial particle velocity distribution is fitted by Maxwellian 
distributions and we find deviations of the high energy tail from the Maxwellian distributions. We estimate the jet's total electron 
content to have a mass between 3.2 x 10^"^ and 1.8 x 10^^ g. Mapping the integrated particle number along the jet trajectory to its 
source region and assuming a typical source region size, we obtain an initial electron density between 8 X 10^ and 5 X 10^^ cm~^ 
that is characteristic for the lower corona or the upper chromosphere. The total kinetic energy of all particles in the jet source region 
amounts from 2.1 X 10^^ to 2.4 x 10^^ erg. 
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1. Introduction 

Polar coronal jets were originally observed by the Extreme- 
Ultraviolet Imaging Telescope (EIT) onboard the So l ar an d 
Heliospheric Observatory (SOHO) (iDomingo et all 1 19951) . 
In white light, the polar jets discovered by LASCO/SOHO 
(1st. Cvr etal.l.ll997h appear narrow and collimated, and expand 
rapidly as they travel through polar regions. They are often asso- 
ciated with an Extreme-Ul traviolet (EUV) jet seen near the solar 
surface (IWang et al.LI 19981) . These jets are often rooted in bright, 
low-lying loop features and are siniilar in appearance to Soft X- 
ray (SXR) jets (IShibata et al.Lll992l:lMoses et al.Lll997h . The es- 
sential acceleration mechanism for all these jets is very likely 
provided by magnetic reconnection. The difference between the 
above and other jet-like features, e.g. chromospheric jets, is the 
altitude where the magnetic reconnection is assumed to occur. 
The higher energy jets tend to be accelerated at a higher alti- 
tude than the lower energy jets (IShibata et al.Ll2007b . The EUV 
and SXR jets are often caused by the reconnection in the upper 
chromosphere or the lower corona. 

Wang et al. (1998) analyzed 27 jets in EIT and LASCO data 
and characterized their motion by three different velocities: the 
leading-edge velocity viead, the centroid velocity Vcen and the ini- 
tial velocity Vinu of the centroid. In all cases Vcen was much less 
than viead, indicating that the jets stretched out rapidly as they 
propagated through the corona. The authors also found that the 
bulk of the jet material decelerated as it propagated from the limb 
to the C2 field of view (EOV). This deceleration was attributed 
to solar gravity. However, the combined results of Vinit < v escape 
and viead > ^escape, and the lack of evidence for downflow in EIT 
and C2 led Wang et al. to propose some in situ acceleration that 



prevents the bulk of the jet material from falling back onto the 
Su n after it was ejected . 

I Wood et al .1 (1 1 999l) improved the jet velocity estimates by 
IWang et al.l (1 19981) using on height- time plots. The authors de- 
termined the trace of the jet centroid from the EIT to the C2 EOV 
and found that the observed kinematic trajectories could be fitted 
with some success by ballistic orbits. They concluded, however, 
that gravity alone was not the only force controlling the jet prop- 
agation. Because o f the similar behavior o f the jets studied, both 
IWood et al.l (Il999l) and lWang et aP (Il998l) suggested that by the 
time jets reached the C2 EOV, they were incorporated into the 
ambient solar wincL 

More recently JKo et al.l (l2005l) studied a jet observed jointly 
by several instruments above the limb. These authors found that 
a ballistic model could explain most of the dynamical proper- 
ties of this jet. In their model it was assumed that the gas was 
ejected upward from the surface with a range of initial speeds. 
The smooth change of the upflow-to-downflow speed at a cer- 
tain altitude derived from the ballistic model was found to be 
consistent with the change of the line intensities from Doppler 
dimming observed by UVCS/SOHO. Owing to a lack of high ca- 
dence coronagraph observations Ko et al.'s study was essentially 
confined to heights below 1.64 r©. 

On June 7, 2007 a big eruptive jet was observed by EUVI, 
CORl, and C0R2 on board the STEREO mission with higher 
spatial and temporal resolution compared to the data from EIT 
and LASCO C2. It extende d from the solar surface to 5 r©. 
The event was also studied by lPatsourakos et aP (l2008l) from the 
stereoscopic viewpoint. They estimated the jet positions and the 
speed of the leading front at different times in the EUVI EOV 
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In this paper we will attempt to analyze the same jet based 
on white-light coronagraph observations at heights beyond about 
1.5 r©. We extend the ballistic approach by Ko et al. ( 2005) by 
quantitatively comparing the density variation from Thomson 
scattered white-light brightness at different heights to the vari- 
ation expected from a ballistic model of the jet particles. This 
method has the advantage that it avoids the estimates of jet cen- 
troids and fronts. These fronts are not well-defined for a diff'u- 
sively spreading plasma cloud, the jet centroids are often difficult 
to determine because a substantial part of jet material is hidden 
behind the occulter. Another advantage of our particle kinetic 
analysis lies is that it also provides a test of the validity of the bal- 
listic model. We therefore receive much more information than 
the conventional leading edge, centroid velocity measurements. 
After the description of the observations in §2, we will introduce 
the ballistic model in §3. In §4 we present the results and try to 
extrapolate our findings to the jet source region and also discuss 
the limits of our model. Finally, we summarize our conclusions 
in the last section. 

2. The data 

2.1. The polar jet in EUVI, C0R1, and C0R2 images 

The jet we investigated was observed over a radial range 
from the solar surface out to five solar radii. We used ob- 
servations from the EUVI and two white light coronagraphs 
CQRl and CQ R2). They belong to the SECCHI instrument suit 
Howard et all 12008). EUVI is a full disk imager with a FOV of 
1.7 r©. CORl and C0R2 are two traditional Lyot coronagraphs, 
with the FOV in the range of 1 .5 r© to 4 r©, and 2.5 r© to 15 r©, re- 
spectively. The angular resolution of one pixel in EUVI, CORl, 
and C0R2 is 1.6 arcsec, 7.5 arcsec, and 14 arcsec, respectively. 

In EUV, the jet could be observed at all four wavelengths 
from the EUVI telescope. The time cadence was 10 minutes for 
304 A and 195 A, and 20 minutes for 284 A. For 171 A, the tem- 
poral resolution was as high as 2.5 minutes. In Figs. [T] and [2] we 
show the diff'erence images of the event at four wavelengths of 
EUVI observed by STEREO A. They were created by subtract- 
ing the pre-event images shown in color in the leftmost column 
of Figs. [Hand [21 The following frames displaying the jet at later 
times is running from left to right with time. Nearly simultane- 
ous images at diff'erent wavelengths are stacked vertically. 

The EUVI observations indicate that this jet contains both 
hot and cool material. The jet eruption first clearly appeared in 
hotter lines 195 A and 284 A at around 05:06 UT, later in 171 A 
at around 05:11, and finally in 304 A at around 05:16 UT. The 
jet life time appeared to be shorter in the hotter lines (171, 195 
and 284 A), and was clearly longer in the cooler line (304 A). 
More detailed descriptions of the EUV observations of this jet 
are provided in lPatsourakos et all (l2008l) . 

The jet we are studying here was also traced in the CORl 
and even in the C0R2 field of view. These observations are the 
main object of our analysis. COR has a polarizer measuring po- 
larized brightness in three directions separated by 120°. From 
this the total brightness is obtained by the standard procedure 
secchi-prep. In Fig.[3]we show the total brightness of the jet ob- 
servations as diff'erence images with respect to the closest pre- 
event image. Again the images are arranged with observation 
time from left to right and stacked roughly synchronously for 
diff'erent telescopes. 

In Fig [4] the jet geometry traced in the FOV of EUVI, 
CORl, and C0R2 is shown by the three thick curve segments. 
A smoothing spline is employed to connect the jet orbit over the 




Fig. 1. The jet time series of diff'erence images at 304 A, 195 A 
and 284 A from top to bottom observed by the EUVI instrument 
on board STEREO A. 



5000 



>- 




2000 



X (Mm) 

Fig. 4. Jet geometry in units of Mm trace from EUVI to CORl, 
then to C0R2 (thick curve segments). The thin solid curve is a 
spline fit to the three traced segments. The half circle indicates 
the solar limb. 

entire height range. In the next section the trajectory of jet parti- 
cles is calculated along this thin solid curve in Fig [4] The good 
alignment of the jet segments from the three instruments indi- 
cates that they are part of the same jet. Notice that there was a 
slight change of about 6 degrees between the jet direction close 
to the surface as seen in EUVI and at altitudes beyond 1.5 r© 
seen in CORl. As a summary, the jet extended out to 5 r©. The 
lifetime in EUVI was around several tens of minutes, in CORl 
the event could be observed for about 1 .5 hours and in C0R2 for 
about 3 hours. 

At 171 A which is the preferred EUV wavelength for plume 
observations, we noticed a preexisting plume to the right-hand 
side of the jet. In CORl, the plume was found cospatial with the 
jet. It is not excepti onal that a pluine is a ligned with a jet in a 
projected 2D image (IWilhelm etal.Ll2011h . However, this more 
or less close alignment along the same line of sight (LOS) may 
be a coincidence. In STEREO B images at 171 A we also found a 
preexisting plume slightly to right-hand side of the jet, in CORl 
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Fig. 3. Jet time series of difference total brightness images observed by CORl (upper row) and C0R2 (bottom row) onboard 
STEREO A. The FOV of each frame in the upper row is 1065'' x 2355'', in the bottom row 2744" x 5614". As a reference, the solar 
radius was around 1002" as observed by STEREO A. 



B the plume signal was too faint to conclude that the plume was 
really close to the jet. In view of the small separation of the two 
spacecraft of only 10 degrees during this event, we still hesitate 
to assume that the jet has a physical relation to the preexisting 
plume. By using difference images, the emission of the plume is 
eliminated for the subsequent analysis. 

2.2. The distance-time brightness relation 

Based on the difference images above, we determined the jet 
intensity in the EUVI images at 304 A and the white-light to- 
tal brightness in CORl and C0R2 images for successive times 
and different distances along the jet axis. The total brightness 
was integrated across the jet cross section for each exposure 
and for each position along the jet axis after the pre-event back- 
ground was subtracted. By integrating the brightness across the 
jet width, at each distance s the density decrease caused by the 
magnetic field line divergence was removed. Note that the inte- 
gration along the line-of-sight direction is implicit in the obser- 
vation of the optically thin jet plasma. The distance s along the 
jet axis remains the only relevant spatial coordinate. The result- 
ing distance-time (DT) total brightness plot is shown in Fig. [5] 
For a clearer view, the brightnesses of CORl and C0R2 in Fig. [5] 
are multiplied by 10^ and 10^^, respectively, to match the inten- 
sity at 304 A. Therefore the absolute values in Fig. [5] are not 
comparable among different instruments. The distance on the or- 
dinate is the length along the jet axis from its footpoint at 1 r©. 
The sampling ds in the direction along the jet axis corresponds to 
the size of one pixel in the respective original images. The time 
is in units of minutes and starts aX = 04:46 UT, which was 
the earliest observation time of the jet in all frames in Figs. [HO 



and [3] It ends at 09:21 UT when we were barely able to see a jet 
signal any more. 

The white-light signal in CORl and COR 2 is caused by 
Thomson scattering at free electrons and therefore the total 
brightness observed by CORl and C0R2 is proportional to 
the line-of-sight integrated column density of coronal electrons. 
Since in Fig.[?]we have subtracted the background and integrated 
the coronagraph signal across the jet cross section, the resulting 
image pixel count B(s, t-t '^jds for CORl and C0R2 in Fig.[5]is 
proportional to the number of jet particles at a distance s along 
the jet in a height range ds resolved by the image pixel, i.e., 

B(s, t-tQ)dsoc N(s, t-t'^)ds, (1) 

where N{sJ - t Q)ds is the number of jet particles in a range 
ds at a distance s from the foot point of the jet and at a time t. 
B(s, t-t 'o) is brightness of the corresponding pixel scaled for the 
different pixel sizes of CORl and C0R2. 

3. Ballistic model 

For the analysis below, we assumed that in the investigated jet 
a package of particles ejected upwards with different speeds vq 
simultaneously during the initiation process. Note that parame- 
ter vo is the velocity component parallel to the jet axis. We did 
not consider the details of the perpendicular velocity, which is 
related to the gyro motion of particles. Moreover, all particles 
that move upward follow the jet axis traced in the images from 
STEREO A. Because we observed the jet very close to the so- 
lar limb, we assume that the projection effect does not effect the 
results much. 



3 



L. Feng et aL: Particle kinetic analysis of a polar jet from SECCHI COR data 



2500 



2000 



J 1 500 



1 000 



500 




50 100 150 200 250 
t - to (min) 



Fig. 5. Image intensity for EUVI 304 and brightness for CORl 
and C0R2 B(s, t-t '^) as a function of time and distance along the 
jet axis, t'^ =04:46 UT, which is a reference time. The intensities 
and brightnesses were integrated across the jet width from the 
difference images (Figs[T] and [3]). The color code on the right 
indicates the 304 intensity in units of number of photons. From 
bottom to top the yellow dotted lines indicate particle trajectories 
with different initial velocities of 250, 300, 350, 400, 425, 450, 
475, 500, 515, 530, 545, 560, 575, 600, 650, and 700 km s'^ 



3.1. The particle trajectory along the jet 

The motion of a charged particle along the unperturbed magnetic 
field is obtained after averaging the forces over the particle gyro 
phase 



ds 
dt 



dv 
dt 



r{sY 



ds 



B(^) - acoii (2) 



s(0,vo) = r© , s(0,vo) = vq. 



Here, s(t - to, vq) denotes the distance along the field line and 
v(t - to, vo) the corresponding velocity in which is the jet ini- 
tiation time. The first term on the right-hand side is the gravity 
force at radius r(s) from the Sun center, is the gravity on 
the solar surface and a is the inclination of the local magnetic 
field that measures the angle between the local radial direction 
and the tangent of the jet axis. The second term represents the 
mirror force driven by the particle's invariant magnetic moment 
jj = mv\l2B and the final term accounts for the deceleration of 
the jet particle by collisions with the background plasma. In the 
following we ignore the mirror and collision term because they 
are small at heights above 1 . 5 r© where the coronagraph obser- 
vations were made. In Sect. 14.31 we justify this choice in more 
detail and discuss possible modifications from our results con- 
cerning the jet properties closer to the surface. 

We solved the second-order ordinary differential equation 
(Equ. [^numerically by a fourth-order Runge-Kutta method. 



Each resulting orbit depends on two parameters, the initial ve- 
locity Vo and the time to that a particle is ejected from the surface. 
Because we assumed that the particles were launched in a unique 
jet event, was chosen to be identical for all orbits. Examples of 
the calculated trajectories for different initial velocities are plot- 
ted as dots in Fig. [5] with vq in the range from 250 to 700 km s"^ . 
The orbit whose apex is close to the edge of the CORl occulter 
at a distance of 500 Mm from the surface has the initial velocity 
400kms-i. 

We have overplotted some particle trajectories in Fig. [5] for 
a somewhat arbitrarily chosen starting time ^o- Qualitatively, we 
find that the brightness of a pixel in the DT-diagram is roughly 
proportional to the number of orbits across it. Because the 
brightness in white-light coronagraphs is caused by the Thomson 
scattering and is hence proportional to the electron density, the 
variation in the jet brightness will be controlled by the particle 
motion. The distance between particles starting at the same time 
with different initial velocities will grow continuously accord- 
ing to their orbit s{t - to, vq) obtained from Equ. [21 To the extent 
that the particle position will disperse, the observed brightness 
of the jet will decrease. This scenario qualitatively agrees with 
the observations in Fig. [5] In the next section, we quantify this 
phenomenon. 

3.2. TheJacobian 

To quantitatively compare the brightness variation caused by the 
ballistic particle motion to the observations from white light 
coronagraphs CORl and C0R2, we estimated the jet particle 
density (Equ.[T]) from the particle trajectories (Equ. [2). For this 
purpose, we equated the number of jet particles N(s, ^ - ^q) in 
distance range s to s-\-ds Sit time t to the number of particles that 
at time to were just in the right range of initial velocities vq to 
Vo + dvo to reach the height range s to s -\- ds at time t. 



N(s,t-tQ)ds = f(vo)dvo 



(3) 



Here, /(vq) is the initial velocity distribution and s = s(t-to, vq) 
depends on the initial particle velocity. Combining Equ.[T]and[3l 
we arrive at 

/(vo) 



B(s,t-tQ)ocN(s,t-tQ) = 



where J(t - to,vo) = 



J(t - to, Vo) ' 
ds(t-to,vo) 
dvo 



(4) 



(5) 



is the Jacobian of the particle orbit. It quantifies the increasing 
dilution of the particle density along the jet with time. In Fig.[5| 
an example of J(t - ^o, vo) is plotted vs time for vo = 500 km s"^ . 
If initially vo deviates slightly, say dvo = 1 km s"\ after a certain 
time, say 50 mins, the diff'erence in distance along the jet is about 
4 Mm because Fig. [6] shows that / at ^ - ^o = 50 min is about 
4 X 10^ s. 

4. Results and discussions 

In this section, the observed brightness distribution B(s, t - t ^) 
along the particle trajectories with given initial velocities vo is 
fitted by the corresponding Jacobians J(t - to, vo) to derive the 
scaling factors /(vo) in Equ.[5]and the optimal jet initiation time 
^0. 



4.1. Fit of the Jacobian to the brightness 

From Equ. 4, we expect the observed B(s, ^ - ^q) to be propor- 
tional to the inverse of the theoretical J(t - to, vo) for every vo if 
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Fig. 7. Comparison of the inverse of the observed brightness B(s, t-t^) times a fitted scaling factor /(vq) (red vertical bars) and the 
theoretical Jacobian J(t - to, vq) (black dashed curve) as a function of time t - after jet initiation at ^o- ^ reference time at 
04:46 UT. The diff'erent diagrams are obtained for diff'erent initial velocities vq, which lead to diff'erent particle orbits s = s(t-tQ,vo) 
along which B(s, t-t '^) is recorded. The vertical range of the red bars indicates the estimated error in B. The width between the two 
solid black curves on both sides of the dased curve shows the variation of / due to an uncertainty in ^o- For more details see text. 
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Fig. 6. One example of the Jacobian as a function of time for 
vo = 500 kms"^ and to = = 04:46 UT. to is the jet initiation 
time and ^ ' is the reference time at 04:46 UT. 



the correct orbit s = s(t - t^, vq) is used in the first argument of 
B and the initial time to is chosen correctly. The constant of pro- 
portionality, /(vo), should then yield the unsealed initial velocity 
distribution of the jet particles. Since /(O, vq) = 0, we rather use 
the inverse of Equ. 4 to obtain estimates for to and /(vq) from 
fits of both sides of the equation. In Eig. [71 we show diagrams 
for diff'erent initial velocities vq of the observed l/B times the 
fitted scaling factor / (shown as red crosses) in comparison to 
the theoretical / (black dashed line) for our best estimates of to 
and/(vo). The initial velocities are chosen to be in the range from 
400 to 650 km s-^ 

The uncertainty in J(t - to, vq) owing to a possible error in 
to is indicated by two solid black curves shifted with respect 
to the central dashed curve by the uncertainty in the initiation 
time. The estimated errors of f/B along the respective trajectory 
s(t - t^, Vo) are indicated by the vertical range of the red bars. 
The error estimate was obtained from the noise in the CORl and 
C0R2 diff'erence images. At each position s along the jet, the 
noise was calculated according to the standard deviation cr of 
the brightness along a circle centered at the Sun center with its 
radius reaching the position s. In Eig. [7] only data points with a 
brightness higher than 3cr are included. The uncertainty of f/B 
indicated by the red bar represents the ±3 cr levels of the bright- 
ness B. 

The slope of J(t - to, vq) at ^ ^ is independent of vq. The 
enhanced slope of / with time reflects the decrease of efl'ective 
gravity with the distance from the Sun. Particles with higher ini- 
tial velocities tend to have less dilution of density (smaller slopes 
for /), which is indicated qualitatively in Eig.[5]as well. The plot 
of curve / in Eig. [7] is terminated either at the time when jet par- 
ticles hit the solar surface in the cases of low vq or at the time 
when the particles leave the EOV in Eig. [5] in the cases of high 

According to Equ. [3 there are two parameters that can 
be modified to obtain a close fit between J(t - to,vo) and 
f(vo)/B(s, t - to). One is the jet initiation time to which is the 
same for all diagrams. A change in to corresponds to a horizontal 
shift of the red bars in all diagrams simultaneously. Next, there 
is a scaling factor /(vq) for each diagram individually. This fac- 



tor represents the relative number of particles with this initial 
velocity in the jet source. Since Equ. \5\ still includes a global 
proportionality constant, the factors /(vq) represent only relative 
number densities at this stage of our analysis. All parameters 
were optimized for a best fit. 

Erom the EUV observations we conclude that the jet started 
between 04:46 UT and 05:06 UT. The first clear jet signal was 
seen at 05:06 UT in 284 A and 195 A. There was a small bright- 
ening at 04:56 UT as well. However, the signal was too weak to 
be identified as the jet initiation. Therefore we assume that this 
jet was ejected in the time range from 04:46 UT to 05:06 UT. 
This uncertainty is also represented graphically in Eig. [7] by the 
two solid curves to either side of the central dotted line repre- 
senting the Jacobian. An initiation time = 04:58 UT gives the 
least sum of the chi- squared deviations for all initial velocities 
Vq. We found a three minute uncertainty, which corresponds to 
the range of to producing slx^ enhancement by 5 % above 

The comparison in Eig. [7] implies that the inverse brightness 
follows the Jacobian curves quite closely. Although the fit is not 
perfect, we can say that the ballistic model in general can ex- 
plain the particle kinematic behavior, and hence the brightness 
variation in this jet. 

4.2. The jet source 

As mentioned above, for each vq the fit of J(t-to, vq) to B(s, t-t ^) 
yields a scaling factor /(vq) that is proportional to the total num- 
ber of particles with this initial velocity. Therefore, the distri- 
bution of /(vo) as a function of vq contains information about 
the initial velocity distribution in the jet source region. This ve- 
locity distribution /(vq) is shown in Eig. [8] together with a fitted 
Maxwellian distribution 



/(vo) = C 



m 

27l€k 



exp 



-m(vo - Vbf 



2€k 



(6) 



Here m is the averaged atomic mass in the solar corona and ap- 
proximated as 1.21 Mp, k is the Boltzmann constant, €k is the 
mean kinetic energy of the particles in the source region and vt a 
velocity shift of the Maxwellian. Because /(vq) is derived from 
the fit of Jacobian to the observed brightness, it does not refer 
to the absolute fraction of jet particles with an initial velocity 
Vq. However, since brightness is proportional to electron number 
density, /(vq) is thus proportional to the normalized velocity dis- 
tribution. Therefore, we added a coefficient C in Equ. [6] In the 
left panel of Eig.[8l the parameters Ek and are chosen to obtain 
the best fit to the observed values of /(vq). The fit results in the 
values of €k = 2.6 x 10"^^ erg for the mean kinetic energy and 
Vh = 230 km s"^ for the velocity shift. 

Because orbits with an initial velocity below 400 kms"^ do 
not reach the field of view of CORl and C0R2, we are lacking 
/(vo) in the velocity range below this threshold. Eor this reason, 
our fit is somewhat insensitive in particular to v^. We therefore 
also considered another fit in the right panel where vt is set to 
zero. We then found a higher mean kinetic energy of 5.4 x 10"^^ 
erg. 

In either case, we found an enhanced higher energy tail at 
Vo > 600 kms"^ that cannot be fitted to a Maxwellian. This de- 
viation from a Maxwellian in the jet source distribution may be 
the result of a specific heating and acceleration process of the jet 
particles. 

In addition to the jet initiation time and the initial velocity 
distribution we have determined the absolute mass ejected by 
the jet. Eor a white light coronagraph image, the electron column 
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Fig. 8. Top: the distribution of the scaling factor /(vq) as a func- 
tion of vo in natural logarithm with a fit by the Maxwellian distri- 
bution. Bottom: similar to the top one, the fit by the Maxwellian 
distribution with fixed = km s"^ . 



density scattering into each pixel is proportional to the calibrated 
pixel total brightness. Since we have subtracted the white-light 
background, only the jet particle density remains from this cal- 
culation. Moreover, the brightness in Fig. [5] was also integrated 
across the jet widths. If we additionally integrate along the dis- 
tance s, we obtain the total number of jet particles visible in the 
coronagraphs at any one time. 

The relations between the image brig htness ^ a nd th e 
electron density w ere established by iMinnae 
van de HulstI (ll95Qh and iBiUrngsl (Il966h . For our purposes, we 
modified them to 



B = 



Bo TicrNedl 



l-u/3 



[(l-u)(2C-A siv?x)]^u(2D-B sm^x)l (7) 



where Bq is the physical mean solar brightness (MSB) that is 
used as units in which calibrated CORl and C0R2 observations 
are expressed. A, B, C, andD are known functions of the distance 
of the scattering location from the Sun's center that express the 
dependence of the scattered polarization on the size of the solar 
disk as seen from the scatterer. u accounts for the solar disk limb 
darkening and cr is the diff'erential Thomson scattering cross sec- 
tion. For u, we used the conventional value of 0.56 for the white- 
light spectral range. Finally, A^^ dl is the electron column density 
along the line-of- sight direction through the jet. As an approx- 



imation, we assumed that the jet was lying in the plane of sky, 
so that we have a scattering angle x = 90° for the Thomson 
scattering. 

For the jet observations at ^ = 05:45 UT when it was most 
prominent, the total brightness has been integrated across the 
width, say integrated over dq to derive B(s, t - to), which again 
was integrated over ds along the jet. If we relate this final inte- 
gration to Equ.[71 we derive J Nedldqds. The result corresponds 
to a total number of 1 . 1 x 10^^ jet particles seen above the occul- 
ter at this time instance. 

From our trajectory analysis this number is related to the 
particles in the velocity range of vq from 400 to 600 kms"^. 
Depending on the extrapolation in velocity space presented in 
Fig. [HI we may extend the above particle number estimate to the 
total number of jet particles. The two cases shown there can be 
considered as two extreme cases. The total number of jet parti- 
cles obtained this way will then lie between 1.6x 10^^ (Fig.[8]left 
case) and 8.9 x 10^^ (Fig.[8]right case). It corresponds to the jet 
mass between 3.2 x 10^^ and 1.8 x 10^^ g. Note that these extrap- 
olations must be treated with some care. They rely on the fact 
that the distribution has a strictly Maxwellian core and even the 
respective parameters are estimates only based on the observed 
distribution of the far tail. Also, the resulting particle number 
may seem large compared to previous estimates merely from 
coronagraph observations, e.g. the COR 1 data at 05:45 UT, be- 
cause it includes all particles involved, also those with a low 
initial velocity vq, which are not seen in the coronagraph images 
at all. 

We may use these number estimates to speculate about the 
particle number density in the source region. If we divide the par- 
ticle number by a typical source volume Vgp, we obtain the prob- 
able number density at the height where the jet heating and/or 
acceleration took place. Here, the volume of this jet source re- 
gion is assumed to be the apparent size of the bright point in the 
195 A pre-event image. We estimate this volume to 2 x 10^^ cm^, 
which yields an electron density between 8 x 10^ cm"^ for the 
case that the jet was heated and accelerated, and 5x10^^ cm"^ 
for the case that the jet was only heated. The densities are typ- 
ical for the upper chromosphere or the low corona and we may 
conclude that the jet material was heated and possibly acceler- 
ated in this height region. Based on the values for Ne, €k and Vbp 
above, the total kinetic energy of all jet particles in the source 
region Ek = (l/2)Ne€kVBp can be estimated which was between 
2.1 X 10^^ erg and 2.4 x 10^^ erg. This energy is consistent with 
typical energies for microflares, indeed our estimate lies in the 
higher energy range of the microflare energy spectrum from 10^^ 
to 10^^ erg. Again we caution the reader that these numbers are 
based on the extrapolation of the particle velocity distribution as 
shown in Fig[8l In particular, we implicitly use vq = kms"^ as 
the lower bound of the distribution for jet particles. The uncer- 
tainty of this lower bound has a strong impact on the estimated 
total number of jet particles, less impact on the kinetic energy in 
the source region. 

4.3. Mirror force and collisions 

The fits of /(vo, to f(v{))/B(s, t) in Fig. [7] as functions of time 
t after jet initiation are in general good but not perfect. A good 
deal of this imperfection can be attributed to image noise, espe- 
cially for high values of t, i.e., low observed brightnesses B{s, t). 
This could be considered as evidence that the particle orbit is suf- 
ficiently well described by the action of the gravity force alone. 
In this subsection we discuss the terms that we neglected from 
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the full equation ([2|, the mirror force and collisions, and find out 
the conditions under which they might become important. 

We approximate the Sun's polar magnetic field by a dipole 
field. Then the magnetic field strength above the pole is 



dH _ ri 



(8) 



where Hq is the field strength at the pole. The mirror force on 
an individual particle can then be written as - jddHldr and 
therefore decreases as oc r"^ with distance from the solar surface. 
Note that ii = mv\ /2H is the particle's magnetic moment and an 
adiabatic invariant of motion provided that the spatial variations 
of the magnetic field are smooth. On the other hand, the gravity 
felt by the particle decreases as oc r"^ and hence less rapidly 
than the mirror force. On the surface of the Sun, comparison of 
the mirror force and gravity force shows that the latter dominates 
as long as the local escape velocity (2^0r|/r(^))^/^ well exceeds 
the perpendicular velocity v^. For a thermal speed on the order 
of the escape velocity of 618 km s"\ a temperature of 46 MK is 
required. 

The other eff'ect neglected in our analysis are possible 
collisions of the jet particles with the coronal plasma back- 
ground. In order to estimate this eff'ect, we used as an esti- 
mate of the collisional deceleration acoii the friction coefficient 
of the Fo kker-Planck colli sion term in a kinetic plasma descrip- 
tion(e.g. Jlshimaru I (Il973h ). The deceleration depends on the ve- 
locity V of the jet particle relative to the thermal velocity Vtherm 
of the coronal background, which is assumed at rest. Then for 
singly charged particles. 



e 

acowiv) = -sign(v) — ^ In A G( 



-), 



^therm 



(9) 



where e is the electron charge. Ad the coronal plasma Debye 
length and In A the respective Coulomb logarithm. An im- 
plicit assumption in Equ. [9] is that the background plasma has 
a Maxwellian velocity distribution. The velocity dependence 
in the Eokker-Planck collision term generally expressed by 
Rosenblu th potentials can then be reduced to the Chandrasekhar 
function (iRosenbluth et al. 1. 1 19571) 

Id erf(x) 

G(x) = - — ( ) 

2 dx X 

The velocity of the jet particles is well ahead of Vtherm on a 
large part of their orbit, therefore G is needed for large ar- 
guments for which it decays as ^ l/2x^. Hence faster parti- 
cles feel less collisional deceleration, which eventually coul d 
lead to a runaway of energetic pa rticles (e.g., iDreice 3 ([T959h : 
ISpringmann and Pauldrach I (Il992h ). 

Eor our estimates we assumed a density as in the left panel 
of Eig. [9] adopted from lOuemerais & Lamyl (l2002h . In the right 
panel of Eig. [9] the resulting deceleration ^coii for diff'erent veloci- 
ties and the gravitational acceleration are compared. The neglect 
of collisions in our analysis is justified above about 1.4 r© for 
particles with 400 kms"^ and for even lower heights for faster 
particles. Eor heights below, the coupling of the jet particles to 
the background plasma is very intense and a single particle anal- 
ysis of the jet seems no longer justified. 

The effect these collisions will probably have is that part 
of the jet's momentum will be transferred to the background 
plasma, which in turn is accelerated. Therefore not all particles 
we see above the occulter may be original jet particles, and the 
distribution function we deduced in Eig. [8] may be the result of 
some interaction of the original jet with the background plasma. 




Fig. 9. Top: Electron density assumed for the b ackground 
corona, adapted from lOuemerais & Lamyl (l2002h . Bottom: 
Collisional deceleration of the jet particles for velocities of 400 
to 600 kms"^ in steps of 50 kms"^ (solid, from right to left) 
compared to the gravitational accelaration (dashed). 



The fact that the distribution in Eig.[8]closely fits to a Maxwellian 
at lower velocities may be evidence for these collisional pro- 
cesses at lower heights. 

In the above collision estimate, we have ignored the solar 
wind in the distribution of the coronal background particles. We 
expect that this does not alter our conclusions substantially be- 
cause below about 1 .4 7?©, where collisions matter, the solar wind 
is probably still subsonic. At greater heights, where the solar 
wind speed becomes significant, the collisional coupling of the 
jet particles is scarce. 

However, a more serious point is that the acceleration mech- 
anism that drives the solar wind may also affect the jet particles 
and should then be added as additional term to the right-hand 
side of Equ. O A physical mechanism for this acceleration has, 
however, not been identified yet and any such acceleration term 
would be highly speculative. 

5. Conclusions 

We have followed the evolution of a big eruptive jet event ob- 
served by SECCHI in both EUVI images and the CORl and 
C0R2 white-light coronagraphs. Based on the distance-time 
brightness analysis, we found that a ballistic model for the jet 
particles in general can explain quite well the brightness varia- 
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tion beyond 1.5 r© in the CORl and C0R2 fields of view with 
gravity as the dominant acceleration of the jet particles. 

Additional parameters were derived or at least estimated, 
such as the initiation time, the initial velocity distribution, and 
the number of the jet particles. The derived initiation time is con- 
sistent with the EUVI observations at lower altitudes. The initial 
velocity distribution was fitted by two Maxwellian distributions 
with diff'erent mean kinetic energies. The good agreement with a 
Maxwellian for lower initial velocities may be due to collisions 
at heights below 1.4 r©. At high initial velocities, the distribution 
deviates from the Maxwellian toward a power law tail that may 
be a result of the jet acceleration process. The total jet particle 
number and kinetic energy sum up to about 1.6 to 8.9x10^^ and 
2.1 to 24x10^^ erg, respectively. 

We neglected the eff'ect of the magnetic mirror force and of 
Coulomb collision. As discussed, they might have some eff'ect 
on the kinetics of the jet particles at lower altitudes. Note that 
these two forces counteract each other: while the collisions with 
the background plasma will decelerate the particles, the mirror 
force accelerates them away from the Sun. Especially the correct 
modeling of the Coulomb collisions below 1 .4 r© requires addi- 
tional assumptions, e.g., about the coronal background density 
and its velocity. 

We outlined the basic idea of a new kinetic jet analysis. In 
the future, a more sophisticated kinetic model of the jet may be 
compared to white-light observations. We have shown that this 
comparison allows one to constrain details of the jet which could 
not be derived in previous studies. More work needs to be done 
in these directions. Moreover, more jet samples are required to 
find out to which extent a jet is embedded in the ambient solar 
wind and how the jet interacts with it. 

PROBA-3, which will be launched in a few years, will have a 
coronagraph with a FOV of 1 .04 to 3 r©. It will provide us with a 
broader initial velocity coverage because the lower velocity limit 
depends on the occulter's size . Therefore we will have less un- 
certainty in the initial velocity distribution, the electron density 
in the jet source region, etc. The higher temporal observations 
with more wavelength coverage from AIA/SDO will help deter- 
mine the jet i nitiation time more p recisely. 

Recently, iRaouafi et al.l (l2008l) found that a jet was very of- 
ten succeeded by a plume above the jet launch site. Interestingly, 
in our jet study a plume was visible in both EUVI and 
CQRl before the jet. This phenomenon was also observed by 
iLites et all (1 19991) and oth er white-light observations in LASCO 
C2 (iLlebaria et al.1. l2002h . However, no definite conclusion is 
given concerning the relation between plume and jet. A time se- 
ries of 3D reconstruction of both plume and jet needs to be made 
to find the answer to this question. 
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